Inelastic hard-rods in a periodic potential 
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A simple model of inelastic hard-rods subject to a one-dimensional array of identical wells is 
introduced. The energy loss due to inelastic collisions is balanced by the work supplied by an 
external stochastic heat-bath. We explore the effect of the spatial non uniformity on the steady 
states of the system. The spatial variations of the density, granular temperature and pressure 
induced by the gradient of the external potential are investigated and compared with the analogous 
variations in an elastic system. Finally, we study the clustering process by considering the relaxation 
of the system starting from a uniform homogeneous state. 
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I. INTRODUCTION 

Recently granular gases, i.e. a large number of macro- 
scopic particles colliding with one another and losing a 
little energy at each collision, have been studied employ- 
ing methods of statistical physics ii*2i2i^ The main prac- 
tical motivation to improve the understanding of the dy- 
namical and static properties of these materials stems 
from their relevance in a variety of technological and in- 
dustrial processesi^ Granular materials also represent a 
new paradigm in the context of non equilibrium statis- 
tical mechanics. In fact, they are open systems which 
can reach a non equilibrium stationary state when the 
energy loss, determined by the dissipative interactions, 
is balanced by the work supplied by an external driving 
force. 

A current problem in this field, motivated by a se- 
ries of recent experiments and theoretical work, concerns 
the behavior of driven granular gases subject to spatially 
non uniform conditionsj These con- 
ditions can be realized experimentally by shaking the gas 
in a container divided in a series of identical compart- 
ments connected by holes. The following phenomenology 
is observed: for strong shaking intensity the grains jump 
from one compartment to the other, whereas for weak 
shaking intensity, they tend to break the discrete sym- 
metry of the system and fill preferentially only some com- 
partments. Such a clustering occurs because the grains, 
unlike an ideal gas of molecules, can dissipate heat when 
colliding. Therefore a spontaneous increase of the num- 
ber of particles within a given compartment entails more 
collisions and a consequent higher energy loss. Thus the 
capability for the particles to leave the compartment is 
strongly inhibited leading eventually to clustering. 

In the present paper we consider a simple modeling of 
the multi-compartment experiment, namely a one dimen- 
sional system of inelastic hard rodsiiiSiiSiSS driven by a 
stochastic thermostat and subject to an external time- 
independent periodic potential. As shown in our recent 
paper fSii the one dimensional geometry of the model not 



only minimizes the amount of computational effort, but 
also allows performing analytical guesses concerning the 
relevant observables. 

As a premise, we recall that in the case of elastic hard- 
rods coupled to a stochastic heat bath and in the pres- 
ence of a periodic external potential most of the equilib- 
rium properties can be computed by using of an exact 
functional relation, proved by PercuSfS^ between an ar- 
bitrary external potential and the equilibrium average 
density profile. It is natural to ask how the behavior 
of the elastic system will be affected upon switching on 
the inelasticity of the collisions. The equilibrium state of 
the elastic hard-rod system is replaced by a correspond- 
ing non equilibrium steady state of the inelastic system. 
The comparison between the properties of the two sys- 
tems offers a useful procedure to understand the interplay 
between the inelasticity, the excluded volume effect and 
the external potential. As we shall see below the model 
shows a variety of non equilibrium states according to the 
values of the control parameters, such as the external po- 
tential, the density and the driving force. These states 
range from spatially periodic configurations to non peri- 
odic configurations in agreement with the experimental 
findings. 

The paper is organized as follows. In section II, we 
introduce the model we employed to study the effects 
of compartments in granular gases and we briefly dis- 
cuss the main feature of the dynamics. In section III, we 
compare the steady state properties of the model against 
the corresponding properties of the equilibrium elastic 
system. In section IV, we consider the relaxation and 
clustering properties of the system. Finally, conclusions 
are presented in section V. 



II. A MODEL OF DRIVEN GRANULAR GAS IN 
COMPARTMENTS 

The main issue of the present paper is to study how 
the properties of a granular fluidized system are affected 
by the presence of a static external non uniform poten- 
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tial. We will show that such a confining potential can en- 
hance the tendency toward cluster formation, generating 
inhomogeneities similar to those observed in some real 
granular systems placed in compartments. The model 
studied consists of a set of N impenetrable particles of 
equal masses, m, and lengths, a, moving on a ring un- 
der the influence of the periodic potential V{x). Two 
mechanisms control the energy of the particles: the in- 
teraction with an heat bath and the collisions among the 
particles. The collisions are instantaneous and binary 
and dissipate a fraction (1 — 7'^)/2 of the total kinetic 
energy of the pair (r being the coefficient of restitution). 
The heat bath, instead, supplies an amount of energy 
which prevents the system to come at rest. In this paper 
we employ, a stochastic energy source, namely the ther- 
mostat previously employed in Refs. (14.18.21 23]. The 
coupling of the particles to the heat-bath consists of the 
combination of a velocity dependent frictional force and 
a random kick. We first write the evolution equation for 
the generic particle i in the system when the collisions 
are not taken into account: 

m^ = -m^^-V'{x,)+m (1) 

where, Xi (i — 1, N), 7 and V'{x) indicate the position of 
particle i, the friction coefficient and the spatial deriva- 
tive of the external static potential, respectively. The 
kicks, (,i{t), are distributed according to a Gaussian law 
characterized by 

and 

The temperature Tb determines the intensity of the heat 
bath. 

The external potential has a spatial period L/A4, 
where L is the system size and M is the number of po- 
tential wells (i.e. number of local minima) in the inter- 
val [— L/2,L/2]. Each well mimics, in our simulation, 
a compartment in Lohse's experiments.^ The boundary 
conditions are assumed to be periodic. We found that 
a rather convenient potential to use in numerical experi- 
ment is that of Refs. 24 25]: 

^ ^ exp{Q;cos(27rxM/L)} 
loia) 

where the zeroth-order Modified Bessel Function, Iq, 
in the denominator, enforces the normalization condition 
1/L dxV{x) = 1. The shape of the potential V{x) de- 
pends upon the parameter a (Fig.^: as it increases, the 
shape of the barriers become more square-like, approxi- 
mating the condition of walls. 

Collisions redistribute the energy and the momenta 
among the particles. For a nearest neighbor pair a col- 
lision event takes place when the separation Xi-^-i — Xi 



equals their length, a and determines a change of the 
velocities of the two particles according to the rule: 

v-^v,-^-^{v^-Vj). (3) 

where the prime indicates post-collisional variables. No- 
tice that in the elastic case, r = 1, the collisions do 
not change the total energy and that the average kinetic 
energy per particle attains asymptotically the value Tb. 
This state of affairs is altered for values of the coefficient 
of restitution less than one and the average kinetic en- 
ergy per particle reaches (for V{x) = 0) a value smaller 
than Tb- 

In the following we will assume a unitary mass to = 1. 
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FIG. 1: Behavior of V{x) for a = 1,2,3, as a increases, the 
barriers approximate the condition of walls and the wells be- 
come more square-like. 

The basic phenomenology of the inelastic model can 
be understood by considering the two competing effects 
controlling the dynamics: the hard core interaction and 
the collisional dissipation. The former favors the equidis- 
tribution of the particles in the various wells, whereas the 
latter promotes the formation of clusters in some wells, 
thus breaking the periodic symmetry of the system. If 
the particles were independent, the passage from one well 
to the next would occur via an activated jump process, 
characterized by a typical exit time, r, proportional to 
exp(Ay/Tf,), where AV — 2 sinh(Q!)//o(a) represents the 
energy barrier between adjacent wells. In other words, 
a local density excess of particles would be washed out 
by a diffusive process whose diffusion constant, D, has 
the Arrhenius dependence, D = Z?o exp(— Ay/Tfc), with 
Do = Tb/m-y. 

In a recent paper— we have shown that, in a granular 
toy model consisting of two wells and only two particles, 
two different jump processes can be clearly identified: 
the free-jump, which occurs when the escape is from a 
singly occupied well, and the correlated-jump that occurs 
when the escape involves a doubly occupied well. The 
resulting scenario is that the jumping dynamics remains 
still Arrhenius-like but the jump-rate strongly depends 
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on the occupation of the well from which the jump occurs 
i.e. 



Tk oc exp 



Tk 



(4) 



where k indicates the occupation number {k — 1 or 2). 
In the multi-well and many particle systems, we shall 
discuss hereafter, the phenomenon is enhanced and the 
escape rate for a given well will depend on the number 
of particle inside that well. 

The steady state properties of our model are inves- 
tigated by means of molecular dynamics simulations of 
Eq. QJ, for an ensemble of = 1024 hard particles of 
length cr — 0.2, in a box of size L — 1024. The heat 
bath parameters are 7 — 0.2 and Tf, = 1.0. The po- 
tential constants a — 1 and M ~ 100, imply 100 wells 
(compartments) of width w = L/M = 10.24, separated 
by an energy barrier of height AV^ — 1.856. The com- 
puter implementation of the granular dynamics has been 
realized by a suitable modification of the event-driven al- 
gorithm to take into account the effects of the potential 
and of the heat-bath. Particle positions and velocities 
within two consecutive collisions are updated according 
to a second order discretization scheme for the dynamics 
Eq. lUJ. 

Simulations start from initial conditions corresponding 
to equally spaced non overlapping particles. After a tran- 
sient time to , the system reaches a steady state character- 
ized by robust statistical properties and one measures the 
generic observable, A(t), by performing the time average 



(A) 



1 



over an appropriate interval, t — 



dTA{T) 



(5) 



III. PROPERTIES OF THE STEADY STATE 

We analyze first the reference elastic system (r = 1) for 
which a true equilibrium state exists. Most of its proper- 
ties are exactly known. For instance, the grand-potential 
ri[p(a;)], from which all the equilibrium properties can be 
deduced, is a functional of the density p{x) whose explicit 
form reads 



L/2 



L/2 



dxp{x) 
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p{x) 



1 — rj{x) 



1 



V{x) 



n 



with /i being the bulk chemical potential of the system 
and 



r]{x) 



dsp{s) 



the local packing fraction. Extremising ^[p] with respect 
to p{x) yields the following Euler-Lagrange equation 



In 



1 - -qix) 



P[p - Vix)] 



p{y) 

dy- 



1 - viy) 



(6) 



whose solution yields the equilibrium density profile. No- 
tice that inserting this solution into one obtains the 
value thermodynamic grand potential of the systemSl 



-n / dx- — — 



(7) 



Using Eq. Q in the uniform limit {V{x) = 0), we obtain 
the bulk pressure —tljL = P = Tbp/{1 — p). 

Equation ^ can also be used to test the average pro- 
file obtained from the numerical solution of the dynam- 
ical equations. Since the potential V{x) varies slowly 
over the hard-core diameter a (i.e. L/M » cr), the 
Euler-Lagrange equation can be solved in the local den- 
sity approximation (LDA) 
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p{x) 


1 — (yp{x) 





1 — (Jp{x) 



(8) 



The density profile Plda{x) solving Eq. (jHJ, is shown 
along with the corresponding quantity obtained from the 
numerical simulation, in figure|3i' and an excellent agree- 
ment is found. In the non uniform case, the pressure is 
obtained from the hydrostatic equilibrium equation: 



dP{x 
dx 



+ p{x)V'{x) = 



(9) 



As we shall see later the grand potential functional 
is related to the probability distribution of numbers of 
particles occupying each cell. 

In the elastic system the temperature and the pres- 
sure can be obtained in two equivalent ways: the first is 
through their kinetic definition the second via their ther- 
modynamic definition. On the contrary in the inelastic 
system both quantities can be obtained only via their 
mechanical definitions. Then temperature is computed 
through the kinetic energy 



N 
i=l 



The pressure is evaluated by considering the impulse 
transferred across a surface in the unit of time^"^-^'* and 
is given by the the sum of the ideal gas term Pid ~ Tgp 
and the collisional Pexc'- 



P = pTg{p) 



a 
Ltob 



Me 



(10) 



fe=l 



Here tob is the observation time, the sum runs over the 
Mc collisions occurring during Tob, and Spk — m6vk rep- 
resents the impulse variation in the fc-th collision. As 
shown in Fig. the pressure profile follows the same 
pattern as the density. Let us remark that with the pa- 
rameters here employed the results of LDA are in fair 
agreement with the corresponding results of the simu- 
lations, thus suggesting that even in the inelastic case a 
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treatment of the excluded volume effect as a local density 
effect might be appropriate. 

Let us consider the effect of the inelasticity. In Figs.Eb 
and 13;, obtained with r = 0.8 and r = 0.6, one observes 
two novel features the temperature profiles become non 
uniform and the density ceases to be periodic. In fact, 
the kinetic temperature decreases in the more populated 
regions, due to the higher collision rate. The particles 
being less energetic, when they belong to a cluster, tend 
to remain trapped even more in the wells. 

Notice that the average pressure is lower than the cor- 
responding pressure in the elastic system at the same av- 
erage density, due to the lower value of Tg with respect 
to Tf). It also varies from well to well, in apparent contra- 
diction with the law of hydrostatic equilibrium. In fact, 
the state described by Figs. 03 and|3; is not stationary, 
but slowly evolves with a characteristic rate associated 
to the average exit time from the wells, much larger than 
in the elastic case. 

At this point, it is interesting to discuss the existence 
of a relationship between the quantities p, P, T, in a gran- 
ular gas, i.e. of an "equation of state". In Fig. Owe plot 
parametrically the numerical values of the pressure ver- 
sus the density obtained by profiles in Fig. [5^ for the 
elastic system. We notice that data are in good agree- 
ment with the exact Tonk's formula™ 



p{p) = n 



p 



1 



up 



(11) 



In a recent paper— we generalized expression to 
the case of a granular system with no external potential, 
and found: 



P{P) = Tg{p) 



P- 



ap'- 



\ — ap 



= T,{p) 



1 — ap 



(12) 



The equation H12|l is equivalent to Eq. Ullfl , upon replac- 
ing Tb with Tg{p). The granular temperature is a func- 
tion of the density since its value is determined by the 
balance between the average power dissipated in the col- 
lisions and the the power supplied by the thermal bath. 
Therefore the larger the number of collisions (i.e. the 
larger the density) the lower the granular temperature, 
which reads 



n 



1 



27 1 — (Tp 



(13) 



Formulae and H12(l have been tested in the homoge- 
neous case and shown to provide reasonable fits to the 
numerical dataiSi 

In the inhomogeneous case {V{x) ^ 0) with inelasticity 
(r < 1), we extracted from the numerical simulations the 
density, temperature and pressure profiles and performed 
the parametric plots shown in Figs.2^,0lD. 

One observes that the plots deviate from the theoret- 
ical predictions and are more scattered. The reason for 
the breakdown of formulae H12|) and can be traced 
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FIG. 2: Profiles of density, temperature and pressure in a 
system of A'' = 1024 rods, size L = 1024, diameter a = 0.2, 
M — 100 potential wells and r = 1.0,0.8,0.6 (a,b,c respec- 
tively), the granular temperature varies with the position 
and is anti-correlated with the density profile. 



back to the inaccurate evaluation of the collision rate in 
the presence of the confining potential. The presence of 
the wells, indeed, strongly increases the collision rate at 
the bottom of each well, thus enhancing the effects of 
inelasticity. 
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FIG. 3: Scatter plot of pressure vs. density for the elastic sys- 
tem. Points are obtained by plotting data of Fig.|2K in a para- 
metric plot {p(x), P{x)} which eliminates the x-dependence. 
Solid line indicates the exact result IIH . 



IV. DYNAMICAL PROPERTIES 

In this section we focus our attention on the cluster- 
ing process which occurs spontaneously, starting from a 
uniform configuration, when the temperature of the heat 
bath is comparable with the energy of the potential bar- 
riers. Since the coUisional cooling determines a decrease 
of the local temperature in the regions of high density, 
some compartments, randomly selected by the dynamics, 
may act as germs for the nucleation of a clustering pro- 
cess. Thus after a long transient, the occupation of few 
compartments may grow at the expense of the remaining 
which become less densely populated. 

An immediate and effective description of the cluster- 
ing phenomenon can be achieved by analyzing the statis- 
tics of the occupation of the wells. We study the probabil- 
ity that the z-th well is occupied by rii particles. In the 
ideal gas, (i.e. independent and identically distributed 
particles) the distribution function of N particles into M 
identical cells follows a Poisson law: 

PM = ^!:^^Eiz^ (14) 

where A = N/AI is the mean (also most probable) occu- 
pation value, and the variance is = pN = N/M = X. 
It is instructive, also in the light of the discussion be- 
low, to link the distribution H14|l to the grand potential 
functional. Let us consider the probability of having a 
configuration with ni,n2, ..,njv/ particles in the M com- 
partments: 

M 

n[K}] = n^K) (15) 

s=l 

For the equilibrium system we write 

n[{n,}] =A/-exp(~r![{nJ]/Tfc) (16) 
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FIG. 4: Scatter plot of temperature a) and pressure b) vs. 
density for a system with an inelasticity coefficient r = 0.6. 
Points are obtained from the profile data of Fig. |5J;. Solid 
lines indicate the theoretical estimate for temperature H13|l 
and pressure 11211 derived in Ref. pH for the system without 
potential. We notice the systematic deviation of the measured 
temperature from the theoretical estimate. Its origin is prob- 
ably due to the large density fluctuations in the boundary 
regions separating the wells. 

where A/" is a normalization factor and indicates the 
ideal-gas coarse grained grand potential 

M 

fl[{ns}] = Tb^^[ns{lnns - 1) - ^Us] (17) 

The external potential, whose role is merely to confine 
the particles into the wells, has been eliminated replacing 
the original microscopic density field p{x) with the coarse 
grained variables n.^. Finally, substituting Eq. (flTjl into 
Eq. HSl) we find: 

mni}]=^l[ ^ (18) 

s—1 

which is identical to expression l|15(l . provided A is re- 
placed by e^^ and introducing a suitable normalization 
factor TV. 
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Such an ideal gas formula has to be modified in the case 
of the elastic hard-rod system, because the hard core re- 
pulsion induces correlations among the positions of the 
particles. The excluded volume constraint can be ac- 
counted for in a simple fashion, by replacing the ideal 
gas functional, Q, by the corresponding hard-rod coarse 
grained functional: 



ilhr[{ns}] = Tb 



M 

E 



Jin? 



l-ln(^„-n,)]-^;;^ 



n 

(19) 

where N„i is the maximum number of non overlapping 
rods which can occupy a single well. The extremization 
of with respect to Ug yields the link between the 
chemical potential /.t, N^n and the mean occupation A 



lnA-ln(Ar,„-A) + 



A 



In this case the distribution reads 



(20) 



(21) 



It displays an enhancement and a narrowing of the peak 
centered at A = N/M, because configurations charac- 
terized by highly populated compartments are unfavored 
(see Fig. The narrowing of the peak can also be un- 
derstood as a consequence of the lower compressibility of 
the hard-rod system with respect to the ideal one. 
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FIG. 5: Occupation probability of the wells for the elastic and 
inelastic system at different values of r = 1,0.9,0.8,0.7,0.6, 
M = 64 wells and N = 512 particles. Elastic curve has a 
Poisson-like shape, given by I2H . peaked around the average 
occupation A = N/M = 8. Curves referring to inelastic sys- 
tem deviate from the elastic one. Inset: Entropic indicator h 
as a function of the inelasticity r. 

The presence of inelasticity introduces a different sce- 
nario with respect to the elastic case. Indeed, the en- 
ergy dissipation, strengthens the correlations and con- 
trasts the hard-core repulsion. Then, as the inelasticity 
parameter r decreases, P{ni) deviates more and more 
from the elastic hard-rod distribution displaying a much 
larger variance, as shown in Fig. |S1 The shape of P{ni) 



reflects the clustering tendency of the inelastic system, 
because a slower decay of P in the region > A indi- 
cates the formation of clusters in few boxes. At the same 
time the peak of the distribution shifts to lower values of 
Tii, showing that the majority of the compartments are 
nearly empty: we observe many wells containing just a 
small number of particles {rii <IC A) and few highly pop- 
ulated wells, in contrast with the elastic case. 

The clustering phenomenon is well represented by the 
following statistical indicator: 



M 

h — ~ y — log — 

i=l 



(22) 



The usefulness of h in characterizing inhomogeneous par- 
ticle distribution in compartments is illustrated by con- 
sidering the two limiting cases. When all particles pile 
up into a single well, /i — > (minimum value), whereas 
h takes on its maximum value, ln(M), when all the M 
wells have identical populations. In intermediate cases, 
the quantity F = exp(/i) represents a measure of the 
average number of occupied compartments. In the in- 
set of Fig. |31 we report the asymptotic value of ft. as a 
function of r referring to a system of = 512 particles 
and M — 64 compartments. The value of h decreases 
monotonically as the system becomes less elastic, indi- 
cating that the grains tend to group together, reducing 
in this way the fraction of occupied wells. Notice that 
as r — > 1 the entropic indicator is slightly larger than 
the value, hn ~ 4.0957 concerning a Poisson distribu- 
tion with the parameter tuned to the simulation setup 
(A — N/M — 8). This discrepancy is a consequence of 
excluded volume effects that result in the narrowing of 
the distribution P{ni) discussed above (Eq, l|2U)- The 
onset of clustering regimes can be emphasized by the time 
evolution of h and by making the direct comparison be- 
tween the elastic and inelastic cases (see Fig. (SJ. The 
system is initially prepared in a state where the grains 
are non overlapping and uniformly distributed among the 
M — 10 compartments. Their velocities are drawn ac- 
cording to a Maxwell distribution at the temperature Tb 
of the heat-bath. The dynamics drives the system toward 
a new set of configurations. 

Figure shows that h (dots) for elastic particles fluc- 
tuates around a constant value which is very close to 
the equi-populated limit ln(M) indicated by the dashed 
line: as expected, there is no evidence of clustering in 
the elastic system. On the other hand, the clustering 
clearly appears from the curves of h{t) referring to the 
inelastic system with r = 0.8. After a transient regime, 
the indicator h{t) becomes stationary, fluctuating around 
a value which is considerably lower than the one corre- 
sponding to r = 1. This means that, due to clustering, 
only few wells result to be effectively populated. The 
time at which the inelastic system departs from homo- 
geneity depends crucially on the number of particles, as 
illustrated in Fig. at A^ = 64 and A^ = 128. We notice 
that the clustering cannot be considered complete be- 
cause h does not vanish but displays small fluctuations 
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around a small constant value. This implies that clus- 
ters exist as dynamical structures, where particles can 
enter and exit. Otherwise, it may even happen that a 
cluster already formed in a well breaks and reforms into 
another well after a long time. However, figure |B1 shows 
that the observable h(t) is a meaningful indicator, be- 
cause it reaches a statistically steady value even when 
p{x, t) is not fully stationary (as discussed in the previ- 
ous section). We stress that if one performs an average of 
p(x, t) on a time larger than the maximum escape time, 
the profile p{x, t) would appear spatially uniform, while 
h{t) for any time t (after the initial transient) indicates 
that the density is non-uniform. Therefore h has to be 
considered the correct order parameter to characterize 
density inhomogeneities in the system. 




30000 



FIG. 6: Evolution of the entropic parameter h, for a system 
starting from a uniform configuration (equi-populated wells), 
with N = M, N = 128 particles of size a = 0.1, M = 11 
wells of width 15.0 (density = 0.775), parameter a = 8.0 and 
energy barriers 3.486. The dashed reference line represents 
the uniform occupation ln(M), the points correspond to the 
elastic system (r = 1), while the lower four curves refer to dif- 
ferent initializations of velocities from a Maxwell distribution 
at temperature Tb — 2, and inlesticity r = 0.8. Thick hor- 
izontal line indicates the average of the asymptotic /i-values 
over 40 runs. 

We turn, now, our attention to a different process: 
we consider the evolution of an initial clustered con- 



figuration, where all the A'' particles at i = are lo- 
cated inside a central compartment in a system made 
of M = 3,5,7,9,15 wellsf^ We expect, according to 
the level of the heat bath temperature T^, two different 
regimes: i) for large the cluster decays toward the fully 
equidistributed state Ui = N/M and limj^+oo F{t) — M, 
since particles tend to fill neighboring cells and distribute 
all over the system; ii) for low Tt the occupation of the 
compartment remains almost constant. It is useful to de- 
fine a characteristic time ig, as the shortest time required 
by the population of the central well, no{ts), to reach the 
fully equidistributed value A = N/M. Unlike the mean- 
field predictions such a lifetime is a strongly fluctuat- 
ing variable, depending on the initial configurations and 
on the "noise" history. Therefore we computed its aver- 
age T = (ts) over 100 different runs and different number 
of wells. The results are reported in figure [7| where the 
average lifetime r is plotted against the inverse of the 
heat bath temperature. We observe the Arrhenius-like 
behavior— with slopes changing with the total number 
of boxes M, in qualitative agreement with the analyti- 
cal predictions The larger the number of wells, the 
smaller the slope, indicating that the most stable cluster 
is the one with only three wells. We also remark that we 
were not able to observe a sudden death of the cluster as 
predicted by theoretical approaches. 
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FIG. 7: Average lifetime r of a cluster initially placed in the 
central well as a function of the heat-bath temperature Tb 
(7 — 0.2) at different number M = 3,5, 7, 9, 15 of wells, each 
of size 20.0. The average has been performed over a set of 100 
randomly selected initial configurations in a granular system 
with A'' — 128 particles of a = 0.1 and r — 0.7, with a = 8 
and potential barrier AV — 3.486 

Among the theoretical treatments of compartmental- 
ized granular gases, the flux-model^iSiiSiiiiSi enjoys a 
great popularity in view of its simplicity. At the heart of 
this approach is the so called flux- function <I>(n), which 
represents a measure of the number of particles crossing 
a section of the system per unit time. The form of the 
flux-function is admittedly phenomenological and follows 
from the experimental observations. The latter require 
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$(n) to be, at fixed driving intensity, an increasing func- 
tion for moderate values of the occupation number, n, 
and a decreasing function for larger values of n, to ac- 
count for the stronger dissipation at high densities. Such 
a non monotonicity entails the possibility of spatially non 
uniform steady solutions.^ 

For the sake of comparison, we extract from simula- 
tions the flux function for our model, by measuring the 
number density of particles, p, flowing outside a given 
well per unit time as a function of the density of par- 
ticles within that well. The simulations start from an 
initial configuration where all the A'o = 50 particles are 
located in the same well of width w = 25 {po = Nq/w) 
and height AV — 3.486. Once the particles reach the top 
of the barrier, they are not allowed to re-enter (absorb- 
ing boundary conditions). As time goes on, we record the 
density, Pout{t), of particles escaping from the well and 
average it over an ensemble of 1000 randomly selected 
initial configurations. The averaging procedure yields a 
curve {pout{t)) from which we estimate its time derivative 
numerically. Finally, a parametric plot of d{pout{t)) / dt 
versus the density of the remaining particles in the well, 
Po ~ {Pout{t)), gives the fiux function <i>(p) of the model. 
It is instructive to compare the properties of the flux 
relative to three different systems: non interacting, elas- 
tically colliding and inelastically colliding particles, as 
shown in figure |S1 In the low density limit the interac- 
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FIG. 8: Flux versus density for a single potential well (M=l) 
with absorbing boundaries. The curves refer to independent 
(dashed), elastic (full) and inelastic (circles) particles. Sim- 
ulations parameters are: initial occupation number A*" = 50, 
o = 0.1, Q = 8, weU size w = 50 (barrier AV = 3.486), bath 
temperature Tb = 1.0 and coefficient of restitution r = 0.8. 
The theoretical prediction of the flux described in the ap- 
pendix is the dot-dashed plot. The inset reports the same 
curves on a larger scale. 

tions are almost negligible and the flux functions display 
a linear behaviors characterized by a common slope. The 
situation becomes more interesting when the well con- 
tains more particles. The flux for non-interacting parti- 
cles grows linearly with p (dashed curve): the larger p. 



the larger is the probability to observe an escape event. 
While, in the system of elastic hard particles, the hard- 
core potential produces two competing effects: on one 
hand the repulsion makes the escape rate for particles 
near the border of the well larger than the corresponding 
rate for non-interacting particles. On the other hand, the 
same repulsion induces a cage effect entrapping the inner 
particles and tends to reduce the escape rate. However, 
from the dot-dashed curve in the inset of flgure|Hl we see 
that the first effect overcomes the second enhancing the 
flux with respect to independent particles. Finally, in the 
inelastic case, the tendency of the particles to form clus- 
ters and to dissipate kinetic energy greatly decreases the 
exit probability from the well, so that the flux is much 
smaller (almost two orders of magnitude) than the flux 
for the independent and elastic systems. We observe a 
slow increase with p and a local maximum clearly visi- 
ble in dot curve of Fig. |S1 The dot-dashed line, in the 
same figure, represents instead the theoretical attempt 
to explain the behaviour of ^{p) through a Kramers' 
theory which combines the excluded volume effect and 
the inelasticity into the escape-time formula (0J) (see Ap- 
pendix) . The theoretical curve is plotted for all the same 
parameters of the simulations except for a unknown pref- 
actor which is used as a fitting coefficient. 



V. CONCLUSIONS 

In this paper, inspired by a series of recent experi- 
ments, we have introduced and studied a model of ran- 
domly driven inelastic hard rods in the presence of an ex- 
ternal spatially periodic potential. The one dimensional 
model presents some of the salient features of higher di- 
mensional models, such as packing, clustering, velocity 
correlations.^"'^ By comparison with the elastic hard rod 
model, used here as reference system, we found that the 
inelasticity has deep repercussions on the configurational 
properties. The fluctuations in density result increased 
and show the tendency of the system to cluster. More- 
over, the (granular) temperature profile ceases to be spa- 
tially uniform showing minima in correspondence of the 
density maxima. As a measure the density fluctuations 
we have introduced an "entropic" indicator, ft, which pro- 
vides information about the statistics of the partition of 
the grains in the various wells. We also performed nu- 
merical experiments concerning the stability of a cluster 
obtained by initially concentrating all the particles in a 
single well. Finally, we considered within the model the 
shape of the flux function which plays a major role in the 
macroscopic description of compartmentalized systems. 
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APPENDIX: FLUX FUNCTION 

In this appendix we present a theoretical argument to 
explain the behaviour of flux function $ measured in sim- 
ulations. The flux of particles outside of a potential well, 
with size w, height AV^ and populated with a density p, 
can be computed via the formula 



dp 
'dt 



Hp) 



P 



rip) 



(A.l) 



where t{p) is the particle escape time across the bar- 
rier. For the present granular system, this time follows 
the Kramers- Arrhenius law for independent particles, but 
with parameters renormalized by interactions. As al- 
ready mentioned, t{p) is the outcome of two compet- 
ing effects, the excluded volume and the granular dis- 
sipation. The dissipation determines the granular tem- 
perature that can be approximated by the solution of 
Eq. 1)13(1 , which holds only for a inelastic system without 
external potential. However for the potential parameters 
chosen in our simulations, formula can be considered 



a fair approximation of the temperature in a well, because 
wells are rather flat and barriers so steep that particles 
feel the confining effect of the potential only in the neigh- 
borhood of well boundaries (Fig|2Jl. The excluded volume 
interaction, instead, is responsible for a mean reduction 
of the potential barrier Ay. This effect is pictorially il- 
lustrated in figureEl where the leftmost/rightmost parti- 
cles (darker), the only ones allowed to jump, are displaced 
from the well center Xc of an amount paw/2. They ex- 
perience an effective energy-barrier roughly estimated as 



dV{p) = V{xc + w/2) - V{xc + p(Tw/2) (A.2) 



and the renormalized Arrhenius law becomes 



r(p) = -exp{5y(p)/T,(p)} (A.3) 



Simulations indicate a very weak dependence of on p 
for the parameters we employed. Therefore, the prefac- 
tor has been considered as a adjustable time scale. 
Numerical solutions of equation H13|) computed at differ- 
ent densities allow to construct the curve of the flux <i> as 
a function of p, reported in figure |S1 in dot-dashed style. 
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